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Abstract 

Representations of quantum state vectors by complex phase space 
amplitudes, complementing the description of the density operator 
by the Wigner function, have been defined by applying the Weyl- 
Wigner transform to dyadic operators, linear in the state vector and 
anti-linear in a fixed 'window state vector'. Here aspects of this con- 
struction are explored, and a connection is established with Gabor's 
'windowed Fourier transform'. The amplitudes that arise for simple 
quantum states from various choices of window arc presented as illus- 
trations. Generalized Bargmann representations of the state vector 
appear as special cases, associated with Gaussian windows. For every 
choice of window, amplitudes lie in a corresponding linear subspace of 
square-integrable functions on phase space. A generalized Born inter- 
pretation of amplitudes is described, with both the Wigner function 
and a generalized Husimi function appearing as quantities linear in 
an amplitude and anti-linear in its complex conjugate. Schrodinger's 
time-dependent and time-independent equations arc represented on 
phase space amplitudes, and their solutions described in simple cases. 
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1 Introduction 



The phase space formulation of quantum mechanics [H El El SI now plays a 
central role in theoretical quantum optics [5] and quantum tomography [6], 
and it has also become an important vehicle for investigations of fundamental 
questions such as the nature of quantization [71 [8] and the quantum-classical 
interface P UHl IH] . New features continue to be revealed by ongoing inves- 
tigations of the underlying structures [121 [T31 E]. 

The overwhelming majority of the very large number of studies in this field 
have focussed on the Wigner function W, the real-valued quasiprobability 
density on phase space F which is the representative of the density operator 
p under the action of the Weyl- Wigner transform W. However, in the last two 
decades a small number of studies [ISl [HI [IZl CSl [13 1201 [21] have gradually 
made clear how the state vector lip), which is after all a more fundamental 
object than p, can also be depicted in the phase space picture, as a complex- 
valued amplitude on F. The recent work of Smith [2T] in particular shows 
that this is achieved most simply by applying W to a suitable multiple of the 
dyadic operator \'ip){ipQ\, where \(po) is an arbitrary normalized state vector 
that, once chosen, is held fixed while {ip) ranges over all states of the given 
quantum system. From this definition it follows that \E' is a multiple of the 
two-state Wigner function [22] M/Vvo' which is thereby given an important 
role in the phase space formulation. 

Several questions concerning this important extension of the phase space 
formulation suggest themselves: 

• In what sense is \4'){'^o\, and hence its image equivalent to \ip) in 
quantum mechanical calculations? 

• What relationships do the Wigner function W and the Husimi function 
[23] have with the phase space amplitude \Ef? 

• What form does the amplitude \Ef take for special system states {ip) 
such as coherent states, or eigenstates of position or momentum, for 
different choices of Iv^o)? 

• For a given \ip), how does the choice of |v?o) influence the structure of 
\E' and its place in the set of functions on F? 

• What forms do the time-dependent and time-independent Schrodinger 
equations take on phase space amplitudes, and how do their solutions 
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look in simple cases, for various choices of Iv^o)? 

• What is the nature of the mapping from \ip) to \E', when viewed as 
a transform of the coordinate space wavefunction ip, and how does it 
relate to other, well-known transforms in the literature? 



These are the questions that we attempt to address in what follows. With 
regard to the last question, we shall see that the transform in question is 
closely related to Gabor's 'windowed Fourier transform' [21], widely used in 
the signal processing literature |25l [26] and also called there the 'short-time 
Fourier transform'. Accordingly we shall refer to the fixed state \^po) that 
appears in the definition of phase space amplitudes \l/ as the window vector, 
window state, or simply 'the window' in what follows. Other names for closely 
related objects in the literature are 'probe functions' [T3] and 'drone states' or 
'fiducial states' [2T]. We shall see also that the choice of an oscillator ground 
state as window gives a that is, up to a factor independent of the state 
the well-known Bargmann wavefunction [27], which is thereby seen as a 
precursor to the more recent efforts [TH [161 HZl [HI [HI |20l [21] to represent 
the state vector in phase space. 

For simplicity of presentation in what follows, we consider mainly one linear 
degree of freedom, ignore spin degrees of freedom, and treat all variables 
as dimensionless, setting Planck's constant h equal to 1. We use hats to 
label operators on the usual complex Hilbert space 1-L of quantum mechanics. 
Variables without hats are defined on F, unless otherwise specified. 



2 Definition of phase space amplitudes 



Our starting point is the observation that pure state vectors can be re- 
placed in all quantum mechanical calculations, without any loss of generality, 
by the corresponding dyadic operators |'?/')((/?o|) for any fixed vector |(^o) of 
unit length, ((^qIv^o) = 1- For example, we can superpose dyadic operators, 

a\i)i) + P\^2) ^ a|^i)(¥'o| + P\^2){^q\ ; (1) 
we can evolve them in time using Schrodinger's equation, 

idtm)) = H\m) ^ ^^t\m){^o\ = i^i^w)(<^oi ; (2) 

we can use them to calculate corresponding pure-state density operators. 
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(3) 



we can use them to calculate expectation values, 

{A) = imij) = Tt (^{\^){^o\)^ A\^){^o\) ; (4) 

and, most important, we can use them to calculate transition amplitudes, 

(^i|z^2) = Tr((|V^i)(y.o|)MV^2)(¥^o|) , (5) 

and not just transition probabilities \ {ipi\il!2)\'^- Transition probabilities, but 
not transition amplitudes, can readily be determined from the density matrix 
or, equivalently, from the Wigner function. 

It is now a simple matter [21] to combine two ideas — the notion of replacing 
state vectors {ip) by dyadic operators \ip){fo\, and the notion of mapping 
operators into functions on phase space using the Weyl- Wigner transform — 
in order to define phase space amplitudes \E'(g,p). We simply set 

<il{q,p) = ^Wm{<fo\){q,p) = ^Tr(\^){^o\A{q,p)) . (6) 

Here A(g,p) is the Weyl-Wigner-Stratonovich kernel operator |28l |29l [30l 
IHTj 1321 ESI EH ES], given in terms of q, p and corresponding (dimensionless) 
canonical operators q, p, by 

A(g,p) =2e2*(P^-«^)n, (7) 
where 11 is the parity operator, satisfying 

gn = -ng, pn = -np, n^ = i. (8) 

We recall at this point that the star product [361 EH EH] of phase space 
functions is defined through the Weyl- Wigner transform by 

W{A) = A , W{B) = B ^ W{AB) = Ai<B, (9) 
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leading to 

{A^B){q,p) = A{q,p)e'''/^B{q,p), {J =d^dp - dpdg) 

= e^(^.^.'-^.^,')/2 A(g,p)i?(g',p')l(.'y)=fep) • (10) 

The expressions involving differential operators are well-defined if A and B 
are polynomial functions. In more general cases, they define asymptotic 
expansions, 



A{q, p) ^ p) = A{q, p)B{q, p) + \iA{q, p) J B{q, p) + 0(2) 



provided A and B represent observables that are asymptotically regular at 
h = Q [3S]. In (11), 0(2) denotes terms of second order in Planck's constant 
(here implicit). 



A string of results now follows from ^ and (|9 11): 



1. Complex phase space amplitudes: For each choice of window vector 
Iv^o), a (distinct) set of amplitudes ^ is defined by ([6]), corresponding 
to the set of all state vectors Each ^ therefore carries implicitly 
a label ipo, which we suppress. For any choice of |<^o), the amplitudes 
^ are complex- valued in general, and carry the same phases as the 
corresponding state vectors up to a constant phase shift determined 

by |<^o)- 

2. Expectation values: Corresponding to (|4]), expectation values can be 
expressed in the form 

{A) = j ^{q,p)^A{q,p)^^{q,p) dT 

= y"{*(g,p)*^(g,p)}A(g,p) dT 

= J ^(g,p) {A{q,p)^^{q,p)}dT. (12) 
Here we have used the result jlOl HI] 
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TtiAB)^ J A{q,p)^B{q,p)dT = J A{q,p)B{q,p) dT , (13) 
for suitably smooth functions A and B. 



3. Transition probabilities: Corresponding to ([S]), transition amplitudes 
and subsequently transition probabilities can be calculated from 



(^i|^2) = J ^iiq,p)MQ,P)dr. (14) 

In particular, 

(V^l^) = Jn^)^{q,p)dT = l (15) 

which suggests that |\E'(g,p)p, like W{q,p), can be regarded as a quasi- 
probability distribution over F in its own right. We shall see below that 
in fact |\E'(g/2,p/2)/2p is a generalized Husimi distribution. 

4. Generalized Born interpretation: Corresponding to (|3]), we have 



W{q,p) = ^{q,p)^<I/{q,p), (16) 



which, bearing ( 12 ) in mind, we refer to as the Born interpretation of \E', 



by analogy with the relation between the wavefunction and probability 



density in configuration space. Then (11) suggests the expansion 



Wiq,p) = ^iq,p)^iq,p) + |z^(g,p) J^(g,p) + . . . , (17) 

which in turn suggests a non-negative approximation to the Wigner 
function, 



W{q,p)^\^{q,p)\^>0. (18) 
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However, (17) and (18) must be treated with caution because \l/ and 



W may not be asymptotically regular at h = [391 B2] . Comparisons 
of |\E'p with W in examples that follow illustrate the difficulty — see 



for example (73) and (74). 



5. Superposition property: For each choice of {ipo), the mapping from 
to \E' is linear and phase space amplitudes can be superposed, 



1^12) = Cil^l) + C2\^2) ^ ^12(g,p) = Ci^i(g,p) + C2^2(g,p) , (19) 

preserving the phase relations between state vectors. 



Note from (16) that if Wi, W2 and W12 are the Wigner functions cor- 



responding to the states and phase space amplitudes in (19), then 



Wu = (Cl^l + C2^2) * (Cl^l + C2^2) 



+C1O2 ^1 * ^2 + C2C{ *2 * 



(20) 



In contrast, it is possible |13] but not straightforward to express ^12, 
for example, directly in terms of Wi and W2- 



Similarly for the generalized Husimi distribution, we have 



1^12! = (Cl^l + C2^2) (Cl^l + C2^2) 



+010^^1^2 + 0201^2^^1. 



(21) 
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6. Normalization and subspace of amplitudes: The results (14) and (15) 
show that, for whatever choice of Iv^o), the complex amplitude be- 
longs to the Hilbert space of square-integrable functions on phase space, 
and is normalized when is normalized. In fact all such amplitudes 
for a given window state lie in a proper closed subspace S^^ of that 
Hilbert space, characterized by the property that 



^{q,p)* W^, (g, p) = m{q,p) /27r , (22) 

where PFI^q is the Wigner function corresponding to the state This 
follows from the identity 



lV')(¥'o|</'o)(¥'o| = |V^)(</'o| (23) 
under the action of W. Note also the identity 

iy^(g,p)*^(g,p) = ^(g,p)/27r, (24) 
which follows from 

|V^)(V^|V')(¥^o| = |^)(¥^o|. (25) 



It can be seen from (14) that if {|'?/'n)} is a complete orthonormal set 
of vectors in "H, then the corresponding set of phase space amplitudes 
{^n} provides a complete orthonormal set in iS^^ under the scalar prod- 
uct given by the RHS of (14). Thus iS^p is a Hilbert space in its own 



right, which for each choice of window provides an image of % within 
the space of square-integrable functions on F. It is clear that "H, 
which after all is commonly realized in the coordinate representation as 
the space of square-integrable functions on the real line, has uncount- 
ably many images in the space of square-integrable functions on the 
phase plane. 



7 



3 A class of integral transforms 

In the coordinate representation, formula ^ takes the form 



(26) 



where if) and 999 are the wavefunctions corresponding to and ly^o) respec- 
tively. 

The first of these formulas shows that, apart from a normalization factor, 
\1/ is the 'two-state' or 'two-sided' Wigner function [31 122] that is sometimes 
denoted W^^p^^{q,p). But as the notation \1/ suggests, we consider it now as 
varying with, and determined by \ip)-, with \(po) held fixed once and for all. 
We emphasize in particular that when \il^{t)) evolves in time, \l/(t) evolves 
also with only the LH member of \ilj{t)){(pQ\ changing (see below). 



The second formula (26) shows that the transform from ijj to obtained 
with any fixed choice of the normalized wavefunction ipo, is closely related to 
Gabor's 'windowed Fourier transform' [211 |25l [26] , which is 



^{q,p) = —= I ^{u)w{u - q) e-''"' du , (27) 
V 27r 



with w{x) the window function. Choosing w{x) = (po{—x) then leads to 



^(g/2,p/2) =2e*'?*'/2$(g,p). (28) 

It is easily checked that \1/ and e*''^/^ $ are also related by a 'symplectic 
Fourier transform', 



e^'?p/2<|.(g,p) = _ / ^(g',p')e^(p''?-5'p)rfr'. (29) 
2n J 

Gabor initially chose a simple Gaussian centred on the origin for the window 
function ipo{x). In the present context, this choice leads to Bargmann's 
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representation [27] of the wavefunction as an entire function on F, regarded 
as the complex plane (see below). 



Note also that the second formula (26) shows each transform in the general 
form 



^(?,P) = y r^o(?>P>w)V^(M)rfM, 

T^^{q,p,u) = (^o(2g-n)e2*^'(«-") . (30) 

From ^ we can also write 

^(g,p) = Tr(|V^)(y;o|A(g,p)) = (</.o|A(g,p)|^) = {%,iq,pM) , 

\%,{q,p)) = A{q,p)\^o). (31) 
In a similar way it can be seen that 

\Q^oiq,P)) = e^(^'^-«^")|^o). (32) 



The transform inverse to (26) is obtained by considering W ^ (^E^) in the 



coordinate representation, which leads to 

^{{x + y)/2,p)e'P^^-yUp = (x|^)(<^o| 



1 

2^ 



= '^{x)My)- (33) 

Hence 



V'(x) = I ^((x + y)/2,p) e'P^'-'-y'^ dp (34) 

27r(y9o(y) 



wherever (po{y) 7^ 0, and also 



^{x) = — II ^iix + y)/2,p)My)e''^'-'^dydp. 



(35) 



Using the definition of W ^ in terms of A{q,p) [35j, the result (35) can also 
be written as 



1 

2^ 



^iq,p)\%oiQ,P))dT 



Similarly 



(36) 



2n 



(37) 



4 Generalized Husimi functions 



The first of formulas (26) says that \E'(g,p) is the Fourier transform of 



Similarly, \E'(g,p) is the Fourier transform of 



{31 



(39) 



and so |\E'(g,p)p is the product of these transforms. By the convolution 
theorem for Fourier transforms [l5], it follows that 



\^iq,p)\' = — I My-z)B,{z)e'^ydzdy 



1 

2^ 



X V^oiq + h- h)(po{q - \z) e*™ dz dy 



IT 



X (po{2q -u- ^v)(po{2q - u + ^v) e du dv . 
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(40) 



Setting 



f{r,v)= MO-\v)Md+\v)e^^Ue. 



(41) 



so that 



MO 



\v) MO+\v) 



1 

2^ 



/(r,w) e"*"^ dr 



(42) 



we then have 



*(g/2, p/2) r = ^ / ^/'(m + ^t;)^(M - ^w)/(r, t;)e-''"(''-") e'^^^ rfr rft; , (43) 



showing that |\E'(g/2,p/2)/2p is a (nonnegative) distribution function from 



Cohen's general class [25l|46]. Note from (28) that 



\^>{q/2,p/2)/2\^d{q/2)d{p/2) = \<^{q,p)\^ dqdp . 



(44) 



In the signals literature, |$(5',p)P is known as the spectrogram of the signal 
ip. It can be regarded as a generalized Husimi function [47j (see below). In 
particular, when the window function is the simple Gaussian corresponding 
to the oscillator ground state (see below). 



we get from (41 ) 



(45) 
(46) 



so that 



|^(g/2,p/2)|2d(g/2)rf(p/2[= \^{q,p)\''dqdp\ =Q{q,p)dqdp, (47) 
where Q is the original Husimi function [23] . 



It is remarkable, especially in view of (17) and (18), that \E'(g,p) T*r \E'(g,p) 
equals the Wigner function for the state \ip), while \E'(g/2,p/2) ^(g/2,p/2) 
defines a generalized Husimi function for that same state, whatever the choice 
of window. In this way \l/ ^ir \Ef and \E' \E' define different aspects of the gener- 
alized Born interpretation of ^ (see above). 
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5 Gaussian windows 



Amongst the simplest choices for a window state is a squeezed state, described 
by a Gaussian function in the coordinate representation, 

= iP'/ny/' e^'^(^-^)/^ e-^'^^-^^'/' , X = xw + tkw/P' , (48) 

where /3, xw and kw are real constants. This has momentum space repre- 
sentation (Fourier transform) 

^^(^^j ^ (^l^ p2^Y/i ^-(k-kw)V2l3^-i(k-kw /2)xw 

= (1/^2^) 1/4 g-/32A(A+A)/4 g-(fc+i/32A)V2/32 ^ ^^g^ 

Then 

|<^o(a;)|2 = e-^'("-^'^)% |v?o(A:)|' = e-('=-'=^)'/^', (50) 

which are peaked at x = xw and A; = kw respectively. Note also that ^^0(2;) 
is the coordinate space representation of a solution of 

-^(/3g+^p//3)bo) = ^|^o). (51) 
Before proceding, we recall that 

n;(g-) = g^g^g...^g = g", y^(^™) = ^ (52) 

and that, following Bopp 



W (Ig") (g,p) = A(g,p) ^g" = [g - i^5JM(g,p) ^ [gBi?]M(g,p) , 

W (Ir) (g,p) = ^(g,p) ^P" = [P + ltd,rA{q,p) ^ [pBRrA{q,p) , 

W (g"l) (g,p) = g"^ A(g,p) = [g + iz9JM(g,p) = MM(g,p) , 

>v(r^) {q,p) =p''^A{q,p) = [p- lz5,]M(g,p) ^ MM(g,p). (53) 
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Applying the first two of tliese results to A = \ilj){(po\ with = 1, it follows 



from (51 ) that every phase space amplitude constructed with the Gaussian 



window state (48) must satisfy 



(54) 



and it is easily seen that the general solution of this equation has the form 



^{q,p) = v^e-^"^/2+^(^+^")^/^G4^) , z = V2{/3q - tp//3) , (55) 
where is arbitrary and we have included an extra z-dependent prefactor 



to for later convenience. The form of G^{z) is determined by ip from (48) 



and the second formula in (26), leading to 



G4z) = (/3Vvr)i/^e^'^(^-^)/^ 



(56) 



6 Generalized Bargmann representation 



Because the exponential factor in (55) is independent of ip, this result estab- 



lishes a generalized Bargmann representation [27] , 



i — > G^{z) , 



(57) 



with G^ an analytic function of the complex variable z (no z-dependence), or 
equivalently of p + i(3'^q, in a Hilbert space Tic with weighted scalar product, 



from (14) and (55) 



{G^, , G^J = (2/7r) / G^,(^)G^,(z)e-^"^+^(^+^)(^+^")/^rf2z. (58) 



Here (Pz = dqdp. The creation and annihilation operators (g =F ip)/\/2 
acting on each lip) in T-L, are represented on as 
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(59) 



and therefore, taking into account the exponential factor in (55), they are 
represented on l-ic as 



{q - ip)/\'2 < — > az - rdz - (3" 



{q + ip) /V2 i — > ad^ - rz + /3^X/V2 , 



a={/3' + l)/2/3, r = (/32-l)/2/3, 



Then 



(60) 



q < — > 



[z + d,)/f3V2 



p i — ¥ i(3{z-d,)/V2-i(3'\. 



(61) 



Note that the RHSs of (60 ) and (61 ) do not involve z and so preserve the an- 
alyticity of any upon which they act. They satisfy the canonical commu- 
tation relations and have the appropriate hermiticity properties with respect 
to the scalar product (58). These expressions generalize those of the usual 
Bargmann representation, which corresponds to the case P = 1 and A = 0, 
that is, the case of a Gaussian window (48) with the same scaling as that 
used for q and p, and centered on x = 0. Then a = 1 and r = 0, and the 



expressions on the RHS of (60) reduce to the familiar z and d^, respectively. 
Furthermore, it can be seen that (56) reduces when /3 = 1 and A 
Definition (2.3) of the Bargmann transform in |2j^ 



to 



The intimate connection between the usual Bargmann wavefunction and the 
Wigner function, and with the phase space formulation of quantum mechan- 
ics more generally, has been discussed previously [IHl EHl EI] from various 
other points of view. Furthermore, it has recently been shown [H] in the 
context of the Gabor transform, that only Gaussian windows give rise to 
spaces of analytic functions of the (generalized) Bargmann type. 



14 



7 A test state and its phase space amplitude 

Consider a normalized test wavefunction 

^v^(l + 2v/2) ^ ^ 

i — > 

^{k) = ^ L-''/'~'' + V2ke-'"/') . (62) 

/0F(l + 2v^) 



which has as Wigner function 



+2\/2(4g2 + p2 _ 1) e-V-pV2 _ ^ [(2g - 1) sin(2p [q + l]/3) 



-2p COs(2p [g + l]/3)] e-(4<72-49+V+l)/3 | (gg) 



and which, for the choice (48) as window with xw = (q) and = (p), 



leads to the phase space amplitude 

^(g,p) = |g-[4/32g2+V+4j(/32-l)gp-4/32(A+l)g]/2(/32+l) 
^g[4i(l-/32A)p+/32(A+l)2]/2{/32+l) ^^^^2 _^ ^^|l/2 
+4^(2/32g - 2zp - P^X) e-[4/3^'^^+V+2i(/32-2),p]/(/32+2) 

^g[4/32Ag+2i/32Ap-/32A2]/(/32+2)^j^^2 _^ 2)3/2 | _ j^g^^j 

Here 

iV = i/4/3/7r(l + 2y2) e^'^^^-^)/^ , A = xvk + ^/ciy/zS' , 

xvi. = (g ) = 1/(1 + 2V2) , fcvK = (p ) = 8V2/9 6^/^73(1 + 2^2) . (65) 
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The values of xw and kw in the definition (48) of the window function have 
been equated to the expectation values (g), 
for the state (62), so that \ipo{x) 



(p) of position and momentum 
and |<^o(^)| are localized near these values 



in position and momentum space, respectively, in accordance with (50). 



Fig. 1 shows the real and imaginary parts of as in (64) when /3 = 1. Fig. 2 
shows the Wigner function W of (63), which is not everywhere positive, and 
the distribution l^'p in this case, for as in (64), with /3 = 1. Fig. 3 shows 
l^'P for P = 0.5 and 2. Comparison of the subplots shows the influence 
of the uncertainty principle: when ^ 1, the structure of better 
delineates the p-dependence of W, while when (3"^ 1, it better delineates 
the g-dependence. 

[Figs. 1, 2 and 3 near here.] 



8 Coherent states in phase space 



The formulas (31) and (32) show that and $ can always be regarded as 



generalized coherent states in Perelemov's sense [52]. However, we now ask 
what phase space amplitudes correspond to the familiar coherent state in 
T-L deflned up to a constant phase by 



1 



^^^{q +ip)\ip^) = -j=\i/j^) , fi = xc + ikceC (66) 
and normalization. Here xc and kc are arbitrary real numbers. It follows 



from (53) that each such state has a corresponding phase space amplitude 
\E'^ satisfying 



1 



V2 



and hence having the general form 



(g+ \idp + ip+ \dq)^^{q,p) 



1 



(67) 



m^{q,p)=K^{w) e 



-ww/2 iiw/\/2 



W 



2{q - ip) . 



16 



Here the precise form of K^j^ is determined by the phase and normahzation 
of and the choice of window Iv^o)- Each choice of \ipo) leads also to 
a corresponding subspace S^p^ of square integrable functions on F through 



(22) and (23), and it is remarkable that the phase space amplitudes (68) 



which we may call coherent states in phase space — must, for every choice of 
window state, form an overcomplete set in the corresponding iS<^„ for varying 
II, just as the coherent states \ip^) form an overcomplete set in "H. 



satisfy (54), leading from (68) to 



If the window is a Gaussian as in (48), then the amplitude \1'^ must also 



K^{w) = const, e 



(69) 



where a and r are as in ( 60 ) . In this case we know also that \E'^ must be of 



the form (55), and we find 



V^{z) = const. e""'/2- g(/.-r/3A)./av^ 



Consistency of (68), (69) and (70) is easily checked, leading to 



(70) 



*^(g,p) = const. e-2[/5^'?^+P^+*(/3^-i)5P-/5^(A+M)^-^(/3^A-;.)p]/(/3^+i) ^ ^^^^ 
and hence 



4/3 



-4[/32(g-,?)2+(p-C)2]/(/32+l) 



vr(/32 + l) 

r] = {xw + xc)/2, C = ikw + kc)/2. 



(72) 



with xw, kw as in (48) and xc, fee as in (66). In (72), we have taken the 



normalization condition (15) into account. On the other hand, the Wigner 



function corresponding to the coherent state \ip^) is [22] 



71 



-{q-^C?-(p-kc? 



(73) 
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If we choose the Gaussian window to be centered on the same coordinate 
and momentum values xq-, kc as the Wigner function, with the same choice 
of length scale — that is to say, if we choose X = fi and /3 = 1 — then 



-2(g-a:c)2-2{p-fec)^ 



TT 



(74) 



which is also centered on the same values as the Wigner function. Note 
however that these two distributions (73) and (74) are not equal, even in this 
case when both are positive. 



Comparing (72) with (73) and (74), we see also the effect of choosing a 
window that is not centred on the key features of the Wigner function. 



9 Schrodinger's equation in phase space and 
evolution of amplitudes 

From (|2|, corresponding to an evolving state vector \ip(t)) in Ti, we have for 
a time-dependent amphtude in phase space, Schrodinger's equation in the 
form 

i dt^iq, p, t) = H{q, p) ^ ^(g, p, t) , (75) 

where H = W(if). Supposing H is not explicitly time-dependent, this equa- 
tion integrates to give 

^{q,p,t) = U{t - to)^^{q,p,to) , U{t) = e'^H^(^'P)*] , (76) 
where the star exponential is defined formally by 

^*[A] ^ j^glj = 1 + ^ + ^ ^ A/2\ + A*A* A/3\ + ... . (77) 

When if is a polynomial in q and p, a more explicit form for the time- 
evolution may be available. For example, if H describes a non-relativistic 
particle and has the form 



H{q,p) = '^p' + V{q) 
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with V a polynomial in q, then from (53) we can rewrite (75) as 



idt'^{q,p,t) = H{qBL,PBL)'^iq,P,t) 



(79) 



leading to 



^0) 



An efficient way to solve (75) or (79) explicitly, when this is possible, is 



to solve Schrodinger's equation in the coordinate representation to get the 



wavefunction il){x,t) and then use that in (26) to construct '^{q,p,t). 
For example, consider the initial-value problem for a free particle 



idtip{x,t) = -ld^ilj{x,t)/dx\ V'(x,0) = v^e-^'^'/Vvr^/^ (81) 
where 7 is a positive constant, with solution 

ij{x,t) = ^e-^'^'/'^^'^'/g{t)7r'/\ g{t) = [1 + z^H]'/', t>0. (82) 

Here the branch of the complex square-root is chosen so that 5^(0) = 1. With 



the Gaussian (48) as window, (82) gives the corresponding time-dependent 



phase space amplitude 



v^(g,p,t) = [4/37M/3W + r)]'/'e 



xe 



-2[p^-f^qHp^gitf+iil3^g{t)^-^^)qp-l3h^M-ip^Xp9m/il3^git)^W) 



as the solution of 



idt<i!{q,p,t) = lpBL^'^{q,p,t) 



corresponding to the initial value 



^3) 



vl/(g,p,0) = [4/37M/3^ + 7')]'/'e 



xe 



2[(5^'y^q''+p^+i{(5^-'y^)qp-(5''j^Xq-il3^Xp]/{l3^+j^) 



15) 
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which may be compared with (71 ) in the case 7 = 1, = 0. 
If the eigenvalue problem 



H{q,p)^<f{q,p) = E<il{q,p) (86) 

can be solved to find a complete set of phase space eigenf unctions '^e„{(1,p) 
— some will be generalized eigenf unctions if H has a (partly) continuous 



spectrum — then a more explicit form of solution to ( 75 ) is 

^{q,P,to) = EnCn^EA(l.P)- (§7) 

Here each sum must be extended to include an integral over the continuous 



spectrum, when appropriate. The eigenvalue problem (86) is distinct from 



the '^-genvalue' problem discussed in the context of the (one-state) Wigner 



function j53l [M], which can be expressed in terms of the solution of (86) 



using (16). 



The expressions (87) are the images under W of the coordinate space formulas 



^(x, t) = E„c„ e-'^"(*-*«) Mx) , ^(x, to) = S„c„ 



where Hijjn = Eni/Jn ■ (88) 

Similarly, orthogonality of the coordinate space eigenfunctions and determi- 
nation of the expansion coefficients, 



1p^{x) i^nix) dx = Smn, Cn = J 1p ^ix) (x , to) dx (89) 

have images 

J^miq,p)'^niq,p)dT = 5mn, = j ^ n{q, p) (q, P, to) dV . (90) 

If the ipn form a complete orthonormal set in the coordinate space repre- 
sentation of Hilbert space, so the form a complete orthonormal set in 
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The simple harmonic oscillator with V{q) = \(f provides the simplest illus- 
tration. In this case, (10) shows that (86) becomes 



W+p') + Wdp-pd,) - \{dl + dl)]^{q,p) = E^{q,p) (91) 
which is easily solved using a phase space variant of the boson calculus. Set 

A=[{q+ lidp) + i{p - lid,)]/V2 = lw + d^, 



(92) 



where we have used checks to distinguish phase space operators, and again 
introduced w = V^{q — ip) as in (68|. 

Solving v4\&o = for the phase space 'vacuum state', normalizing it, and then 
setting = A^'^'^o/y/nJ. for n = 1, 2, . . . , we find 



^^{q,p) = ^^C;:(-1)"-"^«;"^F("-)| 



w e 



-wwjl 



(93) 



m=0 



corresponding to the familiar eigenvalue = n + ^, for ri = 0, 1, 2, . . . . In 



(93), F{yS) and its derivatives are determined by the choice of window and 
the normalization of \E'o; which requires 



F{w)F{w) e-'"'"d'w = 1. 



(94) 



The obvious choice in the present context, ^p^ix) as in (48) with /3 = 1, A = 



leads to F = and (93) is then the Bargmann wavefunction for the 

n-th oscillator eigenstate, apart from the exponential factor. But again we 
emphasize that any convenient window function can be chosen, and every 



choice leads through (16) to the same Wigner functions, which in this case 
are [221 



Ty„(g,p) = (-l)"L„(2[g2+/])e- 



-(g2+p2) 



A, 



(95) 



where L„ is the Laguerre polynomial [SS]. Different choices of y^o lead to 
different phase space amplitudes and hence to different distributions \^n\^ 
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but the same Wigner function (95 ). For example, in the case n = 1, the choice 
^ with /3 = 1, A = gives from ([93| 



vl/i(g,p) = (v^)we-'^»/2 

and hence 

|v]>i(g,j9)|2 = (2/7r) (ww) e"'^'" 
whereas the choice of a 'square' window 



(96) 



(97) 



l/\/2a : |x| < a 
: otherwise 



(9^ 



leads to 

F{w) = (I/tt^/^v^) e'^'/2 Uf (a/y2 - w) + erf(a/y2 + w)] (99) 
in ( [93| ) and from there to 

^i{q,p) = 

(l/vr^/^v^) e-'""'/2 ^(^yj-yj) e'^'/^ [erf(a/y2 - ty) + erf(a/y2 + w) 

+4/v^e-("'+'"')/2 sinh(y2aw)} .(100) 



Fig. 4 shows the Wigner function Wi of (95), and I^E'iP for the Gaussian 
window with /3 = 1, A = as in (97). Fig. 5 shows l^&ip for as in (100) 



in the case of a square window with a = 1. 
[Figs. 4 and 5 near here.] 

10 Eigenstates of momentum and position in 
phase space 



In view of (53), a (generalized) eigenstate of momentum in phase space is 
defined by 
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giving 



PbL ^fco P) = ^O^fco P) ' 



(101) 



'^ko{q,p) = Fkoip) e 



-2i{p-ko)g 



(102) 



with Fk^ip) undetermined. The form of F^^ depends on the choice of Iv^o). 
In fact, when the generahzed momentum eigenfunction (plane wave) 



/V2tt 



(103) 



is inserted into ( 26 ) , it is revealed that 



^koiq,p) = VV^ M^P - h) e-^'^-'^^'^ , 



(104) 



where <^o(^) = (^Iv^o) is the function in momentum space that corresponds to 
(= Fourier transform of) '-poix). It can be checked directly using the integral 
form of the star product [361 EH EHl HI], 



{A^B){q,p) 



\ [ A{q + q2,p + P2)B{q + q3,p + p3)e 



2i(g2P3-'?3P2) 



X dT2dT3 



(105) 



that (16) does indeed hold in this case whatever choice is made for |v?o); in 



every case yielding the (singular) Wigner function 



WkM,P) = Sip-ko)/27r. 



(106) 



In a similar way, the generalized phase space position eigenfunction is defined 

by 



Qbl ^xo(g,p) = a;o^xo(g,p) 



(107) 
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and we find 



(108) 



corresponding to Wxo{q,p) = 5{q — Xq)/2ti. Noting that '^^q in (108) is 
>V(|xo)(v2o|)/v^27r, we then see that the second of formulas (26) is simply the 
image under the Weyl-Wigner transform of the formula 



|^)(V3o| = / \u){u\i)){(pQ\du= / i){u)\u){ipQ\du. 



(109) 



Formulas ( 104 ) and ( 108 ) show again the benefit of choosing a window state 



that is centered on a region of interest in phase space. For example if we 
choose a Gaussian window centered on coordinate and momentum values xw 



tion (104) we get 



and fciy as in (48), (49), (50), then for the phase space momentum eigenfunc- 



^kM.P) = (4/7r3/32)i/4e-(2p-fco-fc^)V2/3^ 



^ ^-H2qp-2qkti-2xwP+'-i:wko+xwkw/'i) 



(110) 



If the window state has ky/ = ko and so is centred on the momentum value 
fco of interest, then 



l*fco(9,P)l' = (4/7r3/32)i/2e-4(P-'=o)V/3^ , (HI) 
which is also centred on p = feg for each value of q, and may be compared 



with the Wigner function (106). 



Similarly, for the phase space position eigenstate ( 108 ) with a Gaussian win- 
dow we get 



^xo(g,p) = (4/3V7r=^)i/^e-'''(2«-"o-"'^)'/2 

y^^i{2qp-2xop~2qkw+xokw+^wkw/'2) 

SO that if the window state has xw = xq, then 



(112) 
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|vl/,„(g,p)|^ = (2/3> 



2 / 3Nl/2„-2/32(g--xo)2 



:ii3) 



centered on q = xq for each value of p. Fig. 
for fco = "~2, with the parameter choices = 



6 shows the real part of \E'fcg 
ko, xw = 4 and /3 = 1 in the 



Gaussian window (48). The imaginary part is similar in form. 



[Fig. 6 near here.] 



11 Oscillator states as windows 



An obvious generalization of the Gaussian window state with wavefunction 
(po{x) as in (48), is the n-th excited state of the oscillator, for some non- 
negative integer n, more precisely, the state obtained by applying a suitable 
'creation operator' n times to the Gaussian 'ground state', to give the nor- 
malized wavefunction 



Mx) = (1/v^) [(/3(a; - A) - a,//5)/y2]>o(x) . (114) 

Windows of this type have been considered recently [56] in the context of 
the Gabor transform; just as the Gaussian choice leads to the (generalized) 
Bargmann transform as noted above, so these oscillator windows lead to 
further generalizations of the Bargmann transform and Bargmann represen- 
tation. 



It follows at once from (114) that if \E'(g,p) is the amplitude corresponding to 
and window state Iv^o), then the amplitude corresponding 



a system state 
to the state lib 



and window state is 



^{q,p) = (1/v^) rnqsR - A) + ipBR/f3)/V2r ^(g,p) 



(115) 



with qBR, Pbr as in (53) 



For example, the coherent state if)^ defined by ( |66| ) leads to the amplitude 



; ^ as in (|71|) when the Gaussian window (|48|) is chosen, whereas the choice 
of the n 



1 ('Mexican Hat') excited state in (114) as window leads through 



(115) to the amplitude 
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V2/3[{2q ~xw- xc) + t{2p - kw - kc)] ^^{q,p)/{(3^ + 1) , (116) 
and hence from ([72| to the non-negative expression 



32/3 



with 1], ( as in (72). This expression should be compared with |^I/^(5',p)P as 
in (72), and with W^{q,p) as in (73). 



12 Concluding remarks 



The introduction of phase space amphtudes extends the phase space for- 
mulation of quantum mechanics, and may be considered to complete that 
formulation by providing images of not only density operators but also state 
vectors. 



The degree of arbitrariness in the definition of the amplitude chosen to repre- 
sent a given quantum state vector reflects the freedom to choose any normal- 
ized state as window state, and may seem surprising at first, notwithstanding 
our final remarks in Section 2. Further study is needed to see how to op- 
timize the choice of window for a given quantum system in a given state, 
for purposes of computation and visualization. Selecting a Gaussian window 
with adjustable width and location in coordinate space and momentum space 
is simple and natural, and often leads to analytic expressions for amplitudes 
representing simple quantum states, as we have seen. This has been useful for 
our purpose here, which was to illustrate basic ideas, but it may not always 
be optimal. A natural choice of window for an evolving quantum system with 
a well-defined ground state and excited states might be its ground state, for 
example. 

Some of the analysis of amplitudes corresponding to quantum states at a 
given instant, mirrors studies of signals in the time- frequency domain done 
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over many years using Gabor's windowed Fourier transform [26l HU [56]. A 
feature of the situation in quantum mechanics that is absent in the case of 
signal processing arises from the very different role that the time variable 
plays in the two cases. In the quantum case, the time is not one of the 
phase space variables, and the description of quantum phase space amplitudes 
evolving in time is a feature absent in the case of signals. Given that the 
window state once chosen is fixed, and does not evolve in time, the choice of 
an optimal window is therefore a more complicated problem in the quantum 
case. 

Other questions that seem worthy of further study in terms of phase space 
amplitudes include the formulation of the uncertainty principle, the descrip- 
tion of interference effects, and the description of quantum symmetries. We 
hope to return to some of these questions. 

Acknowledgements: Thanks to D.R. Stevens for a stimulating conversation, 
and to referees for helpful suggestions. 
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Figure 4: The Wigner function Wi of (95) and as in (97), for a Gaussian 
window having (3 = 1 and A = 0. 
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